grouped analysis

pronoun_by_post_num %>%
  group_by(post_num, total_num_com_in_sub) %>% 
  summarise(
    #agg_author = sum(num_authors),
    agg_i      = sum(num_authors*contains_i     )/sum(num_authors),
    agg_me     = sum(num_authors*contains_me    )/sum(num_authors),
    agg_we     = sum(num_authors*contains_we    )/sum(num_authors),
    agg_us     = sum(num_authors*contains_us    )/sum(num_authors),
    agg_they   = sum(num_authors*contains_they  )/sum(num_authors),
    agg_them   = sum(num_authors*contains_them  )/sum(num_authors),
    agg_you    = sum(num_authors*contains_you   )/sum(num_authors)
  ) %>%
ggplot() +
  geom_point(aes(post_num,agg_i,color=total_num_com_in_sub))
Error in summarise_impl(.data, dots) : 
  Evaluation error: object 'num_authors' not found.

Super Grouped results

all_pronouns <- pronoun_by_post_num %>%
#filter(post_num > 0, total_num_com_in_sub > 100), 
mutate(
  num_til_quit = total_num_com_in_sub - post_num
)
library(plotly)
ggplotly(
ggplot(filter(all_pronouns, num_til_quit < 10), aes(post_num,contains_me,group=num_til_quit, color = num_til_quit)) +
  #geom_point()+
  geom_smooth(se = FALSE)
)
`geom_smooth()` using method = 'loess' and formula 'y ~ x'

  ggplotly(
ggplot(filter(all_pronouns,total_num_com_in_sub < 100, num_til_quit > 6, post_num < 100)) +
  #geom_point(aes(post_num,contains_we,color=total_num_com_in_sub))+
  geom_smooth(aes(post_num,contains_i,color=total_num_com_in_sub))
)
`geom_smooth()` using method = 'gam' and formula 'y ~ s(x, bs = "cs")'

  
  ggplotly(
ggplot(filter(all_pronouns,total_num_com_in_sub < 50, post_num < 20)) +
  #geom_point(aes(post_num,contains_we,color=total_num_com_in_sub))+
  geom_point(aes(post_num,contains_we,color=total_num_com_in_sub))
)

  
  ggplotly(
ggplot(filter(all_pronouns,total_num_com_in_sub < 500, post_num < 200)) +
  #geom_point(aes(post_num,contains_we,color=total_num_com_in_sub))+
  geom_point(aes(total_num_com_in_sub,contains_we,z = post_num))
  #geom_smooth(aes(total_num_com_in_sub,contains_we,z = post_num),se = FALSE)
)
Ignoring unknown aesthetics: z

  filter(all_pronouns,total_num_com_in_sub < 500, post_num < 1000) %>%
  plot_ly(x=.$post_num, z=.$contains_i, y=.$total_num_com_in_sub, type="scatter3d", mode="markers", color = .$contains_we)

  
  filter(all_pronouns,total_num_com_in_sub < 500, post_num < 1000, num_til_quit < 300) %>%
  plot_ly(x=.$post_num, z=.$contains_we, y=.$num_til_quit, type="scatter3d", mode="markers", color = .$contains_we)
ggplotly(
ggplot(filter(all_pronouns,total_num_com_in_sub > 15, total_num_com_in_sub < 400,num_til_quit > 5)) +
  geom_smooth(aes(post_num,ZNcontains_we, color = "We")) +
  geom_smooth(aes(post_num,ZNcontains_i, color = "I")) +
geom_smooth(aes(post_num,ZNcontains_us,color = "Us"  ))+
geom_smooth(aes(post_num,ZNcontains_they,color = "They"))+
geom_smooth(aes(post_num,ZNcontains_them,color = "Them"))+
geom_smooth(aes(post_num,ZNcontains_you,color = "You" ))+
geom_smooth(aes(post_num,ZNcontains_me,color = "Me"  ))
  
)
`geom_smooth()` using method = 'gam' and formula 'y ~ s(x, bs = "cs")'
`geom_smooth()` using method = 'gam' and formula 'y ~ s(x, bs = "cs")'
`geom_smooth()` using method = 'gam' and formula 'y ~ s(x, bs = "cs")'
`geom_smooth()` using method = 'gam' and formula 'y ~ s(x, bs = "cs")'
`geom_smooth()` using method = 'gam' and formula 'y ~ s(x, bs = "cs")'
`geom_smooth()` using method = 'gam' and formula 'y ~ s(x, bs = "cs")'
`geom_smooth()` using method = 'gam' and formula 'y ~ s(x, bs = "cs")'
LS0tCnRpdGxlOiAiMDJfY29tbWVudHNfbnVtYmVyX2FuYWx5c2lzIgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKLS0tCiNncm91cGVkIGFuYWx5c2lzCmBgYHtyfQpsaWJyYXJ5KHRpZHl2ZXJzZSkKcHJvbm91bl9ieV9wb3N0X251bSAlPiUKICBncm91cF9ieShwb3N0X251bSwgdG90YWxfbnVtX2NvbV9pbl9zdWIpICU+JSAKICBzdW1tYXJpc2UoCiAgICAjYWdnX2F1dGhvciA9IHN1bShudW1fYXV0aG9ycyksCiAgICBhZ2dfaSAgICAgID0gc3VtKG51bV9hdXRob3JzKmNvbnRhaW5zX2kgICAgICkvc3VtKG51bV9hdXRob3JzKSwKICAgIGFnZ19tZSAgICAgPSBzdW0obnVtX2F1dGhvcnMqY29udGFpbnNfbWUgICAgKS9zdW0obnVtX2F1dGhvcnMpLAogICAgYWdnX3dlICAgICA9IHN1bShudW1fYXV0aG9ycypjb250YWluc193ZSAgICApL3N1bShudW1fYXV0aG9ycyksCiAgICBhZ2dfdXMgICAgID0gc3VtKG51bV9hdXRob3JzKmNvbnRhaW5zX3VzICAgICkvc3VtKG51bV9hdXRob3JzKSwKICAgIGFnZ190aGV5ICAgPSBzdW0obnVtX2F1dGhvcnMqY29udGFpbnNfdGhleSAgKS9zdW0obnVtX2F1dGhvcnMpLAogICAgYWdnX3RoZW0gICA9IHN1bShudW1fYXV0aG9ycypjb250YWluc190aGVtICApL3N1bShudW1fYXV0aG9ycyksCiAgICBhZ2dfeW91ICAgID0gc3VtKG51bV9hdXRob3JzKmNvbnRhaW5zX3lvdSAgICkvc3VtKG51bV9hdXRob3JzKQogICkgJT4lCmdncGxvdCgpICsKICBnZW9tX3BvaW50KGFlcyhwb3N0X251bSxhZ2dfaSxjb2xvcj10b3RhbF9udW1fY29tX2luX3N1YikpCiAgI2ZhY2V0X3dyYXAofnRvcGljX2xldmVsMSkKYGBgCiMjU3VwZXIgR3JvdXBlZCByZXN1bHRzCmBgYHtyfQoKYWxsX3Byb25vdW5zIDwtIHByb25vdW5fYnlfcG9zdF9udW0gJT4lCiNmaWx0ZXIocG9zdF9udW0gPiAwLCB0b3RhbF9udW1fY29tX2luX3N1YiA+IDEwMCksIAptdXRhdGUoCiAgbnVtX3RpbF9xdWl0ID0gdG90YWxfbnVtX2NvbV9pbl9zdWIgLSBwb3N0X251bQopCgpsaWJyYXJ5KHBsb3RseSkKCmdncGxvdGx5KApnZ3Bsb3QoZmlsdGVyKGFsbF9wcm9ub3VucywgbnVtX3RpbF9xdWl0IDwgMTApLCBhZXMocG9zdF9udW0sY29udGFpbnNfbWUsZ3JvdXA9bnVtX3RpbF9xdWl0LCBjb2xvciA9IG51bV90aWxfcXVpdCkpICsKICAjZ2VvbV9wb2ludCgpKwogIGdlb21fc21vb3RoKHNlID0gRkFMU0UpCikKICBnZ3Bsb3RseSgKZ2dwbG90KGZpbHRlcihhbGxfcHJvbm91bnMsdG90YWxfbnVtX2NvbV9pbl9zdWIgPCAxMDAsIG51bV90aWxfcXVpdCA+IDYsIHBvc3RfbnVtIDwgMTAwKSkgKwogICNnZW9tX3BvaW50KGFlcyhwb3N0X251bSxjb250YWluc193ZSxjb2xvcj10b3RhbF9udW1fY29tX2luX3N1YikpKwogIGdlb21fc21vb3RoKGFlcyhwb3N0X251bSxjb250YWluc19pLGNvbG9yPXRvdGFsX251bV9jb21faW5fc3ViKSkKKQogIAogIGdncGxvdGx5KApnZ3Bsb3QoZmlsdGVyKGFsbF9wcm9ub3Vucyx0b3RhbF9udW1fY29tX2luX3N1YiA8IDUwLCBwb3N0X251bSA8IDIwKSkgKwogICNnZW9tX3BvaW50KGFlcyhwb3N0X251bSxjb250YWluc193ZSxjb2xvcj10b3RhbF9udW1fY29tX2luX3N1YikpKwogIGdlb21fcG9pbnQoYWVzKHBvc3RfbnVtLGNvbnRhaW5zX3dlLGNvbG9yPXRvdGFsX251bV9jb21faW5fc3ViKSkKKQogIAogIGdncGxvdGx5KApnZ3Bsb3QoZmlsdGVyKGFsbF9wcm9ub3Vucyx0b3RhbF9udW1fY29tX2luX3N1YiA8IDUwMCwgcG9zdF9udW0gPCAyMDApKSArCiAgI2dlb21fcG9pbnQoYWVzKHBvc3RfbnVtLGNvbnRhaW5zX3dlLGNvbG9yPXRvdGFsX251bV9jb21faW5fc3ViKSkrCiAgZ2VvbV9wb2ludChhZXModG90YWxfbnVtX2NvbV9pbl9zdWIsY29udGFpbnNfd2UseiA9IHBvc3RfbnVtKSkKICAjZ2VvbV9zbW9vdGgoYWVzKHRvdGFsX251bV9jb21faW5fc3ViLGNvbnRhaW5zX3dlLHogPSBwb3N0X251bSksc2UgPSBGQUxTRSkKKQoKICBmaWx0ZXIoYWxsX3Byb25vdW5zLHRvdGFsX251bV9jb21faW5fc3ViIDwgNTAwLCBwb3N0X251bSA8IDEwMDApICU+JQogIHBsb3RfbHkoeD0uJHBvc3RfbnVtLCB6PS4kY29udGFpbnNfaSwgeT0uJHRvdGFsX251bV9jb21faW5fc3ViLCB0eXBlPSJzY2F0dGVyM2QiLCBtb2RlPSJtYXJrZXJzIiwgY29sb3IgPSAuJGNvbnRhaW5zX3dlKQogIAogIGZpbHRlcihhbGxfcHJvbm91bnMsdG90YWxfbnVtX2NvbV9pbl9zdWIgPCA1MDAsIHBvc3RfbnVtIDwgMTAwMCwgbnVtX3RpbF9xdWl0IDwgMzAwKSAlPiUKICBwbG90X2x5KHg9LiRwb3N0X251bSwgej0uJGNvbnRhaW5zX3dlLCB5PS4kbnVtX3RpbF9xdWl0LCB0eXBlPSJzY2F0dGVyM2QiLCBtb2RlPSJtYXJrZXJzIiwgY29sb3IgPSAuJGNvbnRhaW5zX3dlKQpgYGAKCgoKYGBge3J9CgoKc3RhcnRfcG9pbnRzX2FsbF9udW1zIDwtIGFsbF9wcm9ub3VucyAlPiUKICBmaWx0ZXIocG9zdF9udW0gID09IDEpICU+JQogIGdyb3VwX2J5KHRvdGFsX251bV9jb21faW5fc3ViKSAlPiUKICBzdW1tYXJpc2UoCiAgICBvZmZzZXRfaSAgICAgICAgICAgPSBtZWFuKGNvbnRhaW5zX2kgICApLAogICAgb2Zmc2V0X21lICAgICAgICAgID0gbWVhbihjb250YWluc19tZSAgKSwKICAgIG9mZnNldF93ZSAgICAgICAgICA9IG1lYW4oY29udGFpbnNfd2UgICksCiAgICBvZmZzZXRfdXMgICAgICAgICAgPSBtZWFuKGNvbnRhaW5zX3VzICApLAogICAgb2Zmc2V0X3RoZXkgICAgICAgID0gbWVhbihjb250YWluc190aGV5KSwKICAgIG9mZnNldF90aGVtICAgICAgICA9IG1lYW4oY29udGFpbnNfdGhlbSksCiAgICBvZmZzZXRfeW91ICAgICAgICAgPSBtZWFuKGNvbnRhaW5zX3lvdSApLAogICAgb2Zmc2V0X2F2ZXJhZ2VfdXNlID0gbWVhbihjb250YWluc19pK2NvbnRhaW5zX21lK2NvbnRhaW5zX3RoZW0rY29udGFpbnNfd2UrCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY29udGFpbnNfdXMrY29udGFpbnNfdGhleStjb250YWluc195b3UpLzcKICApCgphbGxfcHJvbm91bnMgPC0gYWxsX3Byb25vdW5zICU+JQogIGxlZnRfam9pbihzdGFydF9wb2ludHNfYWxsX251bXMsIGJ5ID0gYyggInRvdGFsX251bV9jb21faW5fc3ViIiA9ICJ0b3RhbF9udW1fY29tX2luX3N1YiIpKSAlPiUKICAKICAgICAgIyAgICFzdWJyZWRkaXQgJWluJSBjKCJjaGVtaWNhbHJlYWN0aW9uZ2lmcyIsICJwZXJzb25hbGZpbmFuY2UiKSkgJT4lCiAgICAgICAgIyAgc3VicmVkZGl0ICVpbiUgYygiVGhlX0RvbmFsZCIsICJMYXRlU3RhZ2VDYXBpdGFsaXNtIikKICAgICAgICAjKSAlPiUKICBhcnJhbmdlKHNvcnQocG9zdF9udW0pKSAlPiUKICBncm91cF9ieSggdG90YWxfbnVtX2NvbV9pbl9zdWIpICU+JQogIG11dGF0ZSggCiAgICBaTmNvbnRhaW5zX2kgICAgPSAxMDAqKGNvbnRhaW5zX2kgICAvKG9mZnNldF9pICAgKS0xKSwKICAgIFpOY29udGFpbnNfbWUgICA9IDEwMCooY29udGFpbnNfbWUgIC8ob2Zmc2V0X21lICApLTEpLAogICAgWk5jb250YWluc193ZSAgID0gMTAwKihjb250YWluc193ZSAgLyhvZmZzZXRfd2UgICktMSksCiAgICBaTmNvbnRhaW5zX3VzICAgPSAxMDAqKGNvbnRhaW5zX3VzICAvKG9mZnNldF91cyAgKS0xKSwKICAgIFpOY29udGFpbnNfdGhleSA9IDEwMCooY29udGFpbnNfdGhleS8ob2Zmc2V0X3RoZXkpLTEpLAogICAgWk5jb250YWluc190aGVtID0gMTAwKihjb250YWluc190aGVtLyhvZmZzZXRfdGhlbSktMSksCiAgICBaTmNvbnRhaW5zX3lvdSAgPSAxMDAqKGNvbnRhaW5zX3lvdSAgLyhvZmZzZXRfeW91ICktMSkKICAgICkgJT4lIHVuZ3JvdXAoKQoKZmlsdGVyKGFsbF9wcm9ub3Vucyx0b3RhbF9udW1fY29tX2luX3N1YiA8IDUwMCwgcG9zdF9udW0gPCAxMDAwLCBudW1fdGlsX3F1aXQgPCAzMDApICU+JQogIHBsb3RfbHkoeD0uJHBvc3RfbnVtLCB6PS4kWk5jb250YWluc193ZSwgeT0uJHRvdGFsX251bV9jb21faW5fc3ViLCB0eXBlPSJzY2F0dGVyM2QiLCBtb2RlPSJtYXJrZXJzIiwgY29sb3IgPSAuJFpOY29udGFpbnNfd2UpCgoKZ2dwbG90bHkoCmdncGxvdChmaWx0ZXIoYWxsX3Byb25vdW5zLHRvdGFsX251bV9jb21faW5fc3ViID4gMTUsIHRvdGFsX251bV9jb21faW5fc3ViIDwgMTAwLG51bV90aWxfcXVpdCA+IDUpKSArCiAgZ2VvbV9zbW9vdGgoYWVzKHBvc3RfbnVtLFpOY29udGFpbnNfd2UsZ3JvdXA9dG90YWxfbnVtX2NvbV9pbl9zdWIsY29sb3IgPSB0b3RhbF9udW1fY29tX2luX3N1Yiksc2UgPSBGQUxTRSkKKQoKCmdncGxvdGx5KApnZ3Bsb3QoZmlsdGVyKGFsbF9wcm9ub3Vucyx0b3RhbF9udW1fY29tX2luX3N1YiA+IDE1LCB0b3RhbF9udW1fY29tX2luX3N1YiA8IDQwMCxudW1fdGlsX3F1aXQgPiA1KSkgKwogIGdlb21fc21vb3RoKGFlcyhwb3N0X251bSxaTmNvbnRhaW5zX3dlLCBjb2xvciA9ICJXZSIpKSArCiAgZ2VvbV9zbW9vdGgoYWVzKHBvc3RfbnVtLFpOY29udGFpbnNfaSwgY29sb3IgPSAiSSIpKSArCmdlb21fc21vb3RoKGFlcyhwb3N0X251bSxaTmNvbnRhaW5zX3VzLGNvbG9yID0gIlVzIiAgKSkrCmdlb21fc21vb3RoKGFlcyhwb3N0X251bSxaTmNvbnRhaW5zX3RoZXksY29sb3IgPSAiVGhleSIpKSsKZ2VvbV9zbW9vdGgoYWVzKHBvc3RfbnVtLFpOY29udGFpbnNfdGhlbSxjb2xvciA9ICJUaGVtIikpKwpnZW9tX3Ntb290aChhZXMocG9zdF9udW0sWk5jb250YWluc195b3UsY29sb3IgPSAiWW91IiApKSsKZ2VvbV9zbW9vdGgoYWVzKHBvc3RfbnVtLFpOY29udGFpbnNfbWUsY29sb3IgPSAiTWUiICApKQoKCgoKCgoKCiAgCikKCmdncGxvdGx5KApnZ3Bsb3QoZmlsdGVyKGFsbF9wcm9ub3Vucyx0b3RhbF9udW1fY29tX2luX3N1YiA+IDE1LCB0b3RhbF9udW1fY29tX2luX3N1YiA8IDQwMCxudW1fdGlsX3F1aXQgPiA1KSkgKwogIGdlb21fc21vb3RoKGFlcyhwb3N0X251bSxjb250YWluc193ZSwgY29sb3IgPSAiV2UiKSkgKwogIGdlb21fc21vb3RoKGFlcyhwb3N0X251bSxjb250YWluc19pLCBjb2xvciA9ICJJIikpIAogIAopCmBgYAoKCgo=